Nonequilibrium charge density wave ordering from anomalous velocity 

in itinerant helical magnets 
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The Karplus-Luttinger anomalous velocity is shown to lead to electric-field-induced charge ac- 
cumulation in nearly ferromagnetic noncoUinear magnets with itinerant electrons, like MnSi. For 
helical magnetic ordering, the balance between this accumulation and the Coulomb interaction leads 
to a nonequilibrium charge density wave state with the period of the helix, even when such accu- 
mulation is forbidden by an approximate gauge- like symmetry in the absence of electric field. We 
compute the strength of such charge accumulation as an example of how unexpected many-electron 
physics is induced by the inclusion of the one-electron Karplus-Luttinger term whenever the local 
exchange field felt by conduction electrons does not satisfy the current-free Maxwell equations. 
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Recently the Karplus-Luttinger "anomalous velocity" 
in crystal structures without inversion symmetry has 
been the subject of renewed interest. It has been reinter- 
preted as an additional term in the semiclassical equa- 
tions of motion resulting from a Berry phase in momen- 
tum space^, and shown to be determined solely by elec- 
tronic properties at the Fermi surface^, consistent with 
the spirit of Landau's Fermi liquid theory. A natural 
question is whether the existence of a new term in the 
single-electron equation of motion modifies standard in- 
teracting electron physics. We answer this question in 
the afRrmative for long-period itinerant helimagnets like 
MnSi and related systems: the anomalous velocity leads 
to a nonequilibrium steady-state charge density wave 
(CDW) transverse to an applied electrical field. 

Such ordering is forbidden by symmetry in the ab- 
sence of an electrical field, but the combination of an ap- 
plied spatially constant electrical field with the spatially 
variable effective "Hund's rule" magnetic field yields a 
charge density wave whose amplitude is balanced by the 
Coulomb interaction. Similar noncoUinear orderings may 
exist in other systems like dilute magnetic semiconduc- 
tors, where simple estimates predict a much stronger 
version of the effect we find. We show that the same 
type of nonequilibrium charge ordering is generic to non- 
coUinear magnets that support effective magnetic fields 
(from exchange interactions) that are inconsistent with 
the current-free Maxwell equations. 

As in the two-dimensional quantized Hall effect, there 
are finite transverse conductances in three-dimensional 
transport that do not depend on relaxation and hence 
are "intrinsic." Here we focus on the Karplus-Luttinger 
anomalous velocity in magnetic materials with broken 
time-reversal symmetry ("T-breaking materials"), which 
has been argued to explain the anomalous Hall effect in 
ferromagnetic semiconductorsi. A similar intrinsic trans- 
verse conductance underlies the "spin Hall effect" pre- 
dicted for some doped semiconductors without broken 
time-reversal symmetrySiii^. Current in all these situ- 
ations flows perpendicular to electric field and another 



vector: in the ordinary Hall effect, the other vector is 
magnetic field; in the anomalous Hall effect, it is the lo- 
cal magnetic moment; and in the proposed intrinsic spin 
Hall effect, it is the spin polarization that determines the 
flow of the current. Because the direction of the current 
is determined by this third-party vector, it is interesting 
to ask what will happen if this vector is not homogenous. 

The standard derivation of the anomalous Hall effect 
contribution from the Karplus-Luttinger term depends 
on nearly constant local moments, so the nearly ferro- 
magnetic noncoUinear magnets like MnSi^t^ are natural 
candidates. An intuitive understanding of the induced 
CDW we flnd in such materials can be gained by con- 
sidering a stacking of several blocks of an itinerant fer- 
romagnet (e.g., iron), in such a way that the local fer- 
romagnetic direction rotates in a helix (Fig. 1). Then 
under an applied electric field perpendicular to the he- 
lix direction, the different Hall currents in each block 
will lead to a buildup of charge at the junctions between 
different blocks. The result we find for itinerant heli- 
magnets can be viewed as a continuum version of the 
above physics, in which the transverse current fiow and 
charge buildup is dominated by the anomalous Hall ef- 
fect in materials with broken time-reversal and inversion 
symmetries. This charge buildup will be largely com- 
pensated by screening by the conduction electrons, but 
there is a residual lattice distortion and charge oscilla- 
tion in the helix direction that should be experimentally 
detectable. 

Long-period modulations (periodicity over several 
hundred angstroms) are also found in SrFeO^ and 
CaFeO^. We find that in these helimagnets and some 
other configurations, an external electric field in certain 
directions will induce a charge density wave, and that 
this nonequilibrium charge ordering is directly related to 
the symmetry of the magnetization. The estimates we 
give of the strength of this charge density wave are easily 
carried over to other materials; as the anomalous velocity 
only depends on the electrons or holes close to the Fermi 
surfac^, the detailed band structure is irrelevant. 
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Since we will consider materials with no symmetry un- 
der spatial inversion and with nonzero spin-orbit cou- 
pling, the leading spin-orbit coupling should be linear in 
electron momentum. The free-electron Hamiltonian is 

^ = II4.(ek + 7k-a)ck.. (1) 

For simplicity, we will later approximate the dispersion 
relation by ek = ?i^/(2me) • fc^. The detailed structure of 
the spin-orbit coupling term is not too important for the 
following calculation, and in fact the same effect appears 
for the Luttinger Hamiltonian for a spin-3/2 band with 
only (k • s)^ terms. The modified semiclassical equation 
of motion is 

hk = —eE — ex X B 



The Berry curvature 17 can be regarded as a magnetic 
field in momentum space: in terms of the periodic part 
Un of the wavefunction for Bloch band n, its form is 

aKfc) = .,,Im^|g||^y (3) 

We can easily see in that the helicity is a conserved 
quantity, and for the same momentum, electrons with 
different helicities have different energy, so actually here 
we have two bands, with two different helicities. These 
two bands have a degenerate point at k = 0, close to this 
point, using Q, we will see the Berry phase looks like a 
monopole in momentum space, 

a:(A,fc)-A^. (4) 

Here A is the helicity k • a. Although we use this form 
for the Berry phase in our calculations, the specific form 
of the Berry phase is not important: the anomalous Hall 
effect is basically a Fermi surface effect. If there is a 
gap between two bands, an infrared cut-off at fc = is 
necessary, but this will not contribute to the final result. 

Under an electric field E, the solution of Q has both 
normal acceleration along electric field and an anomalous 
velocity in the plane normal to E, leading to an anoma- 
lous Hall current^ 

J = E^ fd'mi\)xE)fie). (5) 

Here /(e) is the distribution of particles. If /(e) is just 
the equilibrium distribution, the anomalous current is 
zero because the equilibrium distribution is even under 
k — > — k, and the Berry phase is odd. More generally, 
the anomalous Hall effect is forbidden by time-reversal 
symmetry: electric field is even under time reversal, but 
current is odd, and because the current is perpendicular 
to electric field, this anomalous Hall current has no dis- 
sipation, so without another vector which is odd under 



time-reversal there can be no Hall current. In magnetic 
materials the ordered moments provide this additional 
vector, analogous to an external applied field. 

The ordered moments produce the Hall current via the 
"local Hund's rule" coupling to spin, whose origin is in 
exchange interactions rather than magnetic dipolar in- 
teractions. This interaction between itinerant electrons 
or holes and the local magnetic order can be modelled in 
a simple Hamiltonian 

H' ^hm- s. (6) 

Here m is the local magnetic moment. It is assumed 
that the background moments are strongly ordered so 
that it is appropriate to use the expectation value of m 
as the Hund's rule field and ignore fluctuations. We will 
assume for now that the magnetic configuration has con- 
stant magnitude and normalize |mp = 1. 

The conduction electron spin density s includes spin- 
1/2 matrices in MnSi and spin-3/2 for the Luttinger 
model of GaAs (here we define s to be dimensionless 
so that h has units of energy). The exchange field h is 
proportional to the magnitude of the ordered local mo- 
ment. If Mn spins are fully polarized, h ~ NMnSJ<^, 
where Nmu is the density of Mn ions, S is the spin of 
Mn ion, which is 5/2, and J is the exchange energy be- 
tween local moment and itinerant electrons, which we 
take to be 50 meV nm^— . In MnSi, the Mn ions are 
not fully polarized and the local moment is about 0.4/iB 
per ion^^. So in this material h is suppressed to 8% of 
its maximum value. In MnSi, the lattice constant is 4.6 
angstroms and there are 4 Mn ions per unit cell, so we 
obtain Nmu — iO/nm? . Collecting all the results above, 
we finally estimate h = 300 meV in MnSi. 

On symmetry grounds, if both time- reversal and inver- 
sion symmetries are broken, one expects a Hall current 
from the anomalous velocity term: 

Jah = ctahE X m. (7) 

Both sides of this equation are odd under time reversal, so 
a nonzero anomalous Hall effect is consistent with broken 
time-reversal and inversion symmetries. 

In this problem we have three different energy scales: 
the Fermi energy, spin-orbit coupling energy, and the 
Hund's rule coupling between itinerant spin and local 
moment. For convenience we assume Ef 3> X'^kf ^ A/i, 
although our main result does not depend on this as- 
sumption. In this case the basic band structure of this 
material is still applicable and we can regard the effect of 
local moment as a perturbation. Suppose that m is along 
the z direction. Then the Fermi sea with positive helic- 
ity will be pushed down along — z, because if an electron 
with momentum along negative z has positive helicity, 
it will benefit from the interaction with local moment. 
For the same reason, the Fermi sea with negative helicity 
will be pushed up along z. Thus the Fermi surface has 
an anisotropic Fermi wavenumber kp. 
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Because the energy at the Fermi surface must be con- 
stant, we have 



7, 



2m, 



ki + X-ke + X- cost 

2m, " 2 2 



Solving these equations, we obtain 



hmgX cos 9 



h^k 



'f 



(8) 



(9) 



We have used the assumption that the Fermi energy is 
much greater than the spin-orbit couphng to ignore the 
7 in the denominator. So the unknown value of 7 is irrel- 
evant as long as it satisfies the assumption above. This 
is analogous to the situation in the Rashba model, which 
in the absence of interactions and disorder has a uni- 
versal spin-Hall conductivity independent of spin-orbit 
coupling^. Here kj is the Fermi vector at 9 = tt/2, is 
the electron effective mass, 9 is the angle between k and 
m, and A is the helicity. Now because of this anisotropy, 
integrating over the Fermi sea as in and summing over 
the two helicities finally obtains a nonzero Hall current: 
the transverse conductivity is 



e^hme 1 
kf' 



(10) 



For a numerical estimate, take 1/fc/ ~ 1 A and replace 
the effective mass by the bare mass: then the Hall con- 
ductivity a AH — 2 X 10'^ fl~^cm~^ . In very clean samples 
of MnSi, the diagonal conductivity a is about 10 fifl^^ 
cm~^ below 20 K, the temperature below which the local 
moments exis1i&, so (jahI<^ ~ 10~*. As discussed below, 
this ratio will increase in slightly disordered samples. 

This l/fcy behavior appears strange because it predicts 
that as the particle density vanishes, the Hall effect will 
be infinitely large. However, the above calculation is 
based on the assumption that the coupling with local mo- 
ment is much weaker than the Fermi energy. If fc/ is very 
small, the Hund's rule coupling dominates, almost all the 
spins are aligned along the direction of local moment, and 
it is the spin-orbit coupling that is a perturbation. 

All the calculations above are based on the nonvanish- 
ing Berry phase, which depends on adiabaticity of elec- 
tron motion. Heuristically, because At • Ai? > ?i/2, the 
relaxation time Ai should be large in order to make sure 
Ai? is very small compared with the band gap at Fermi 
surface. Although the precise strength of spin-orbit cou- 
pling in the materials we study is unknown, it has to 
be much bigger than ft., the coupling to local moment, 
for the perturbative calculation above to be valid. So if 
T » ft/ft, « 10~^^s, the calculation should be valid: this 
is a typical relaxation time in a metal in room tempera- 
ture, so at low temperature our assumptions are satisfied. 

Now consider what will happen if the local moment is 
not perfectly coUinear. First we discuss the helical case 
which is relevant to MnSi and other materials discussed in 
the introduction^^. For simplicity, we take local moments 




FIG. 1: Geometry of charge accumulation in a helical mag- 
netic background. In an electric field E perpendicular to the 
helix direction, a transverse CDW forms with the same mod- 
ulation vector q as the helix. 



The modulation period 21x1 q is typically at least 100 A, 
which means the local moment is rotating very slowly 
when we move along x. Because it is rotating very slowly, 
particles in different regions of the material only see a 
homogenous local moment, and locally have the same 
momentum distribution as previous several paragraphs. 
First keep the electric field along y: then based on {T)) 
and H10|l . the local Hall current density is 



JahIx) = cos{qx)aAHEy 



(12) 



However, electric charge is always a conserved quantity, 
and obviously, this anomalous charge current does not 
satisfy the static continuity condition V ■ J ah = 0, so 
charge will accumulate, inducing an electric field. Finally, 
given some dissipation mechanism, the system will reach 
a steady state in which the diagonal current from the 
induced electric field cancels the original Hall current: 



aE' = -Jah- 



(13) 



Here a is the normal diagonal conductivity of the mate- 
rial. From the Maxwell equation V • D = Anp, we obtain 
the relation between accumulation and electric field 



p{x) — gsin(gx) ('i') AH 



(14) 



xn(x) = cos{qx)z + sm{qx)y. 



(11) 



e(q) is the static dielectric function at wave vector q. This 
p{x) describes the bare or external electric charge den- 
sity, which will be strongly screened if e{q) ^ 1. 

Let us estimate the magnitude of this charge density 
wave effect. First, in metal e(g) is given by e(g) = 
1 + Qtf/i^^ if = 10 get e(g) ~ 2 x 10^. Also 

we have aAn/o' ~ 10^'*. If we take electric field to be 
100 V cm~^, the density modulation is about lO^^cm"'^, 
which because of its periodicity should be measureable 
in a diffraction experiment. Furthermore, from H14|) the 
density modulation is proportional to the ratio between 
anomalous Hall conductivity and diagonal conductivity. 
Reducing the diagonal conductivity by adding impuri- 
ties will lead to increased charge modulation until the 
relaxation time becomes so short that the Berry's phase 
assumption discussed above is violated. 

This charge density wave exists for helical direction x 
when the electric field has in any component in the yz 
plane. However, if the field is along x, the Hall current 
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will be in yz plane, and it is homogenous in this plane, 
so V • J = and there is no accumulation. Actually, 
the formation of charge density wave is because of the 
breaking of an approximate spiral symmetry possessed 
by the long-period helical spin state. A homogeneous 
background has 3 translational symmetries, Tx, Ty and 
Tz ■ When the helical magnetic moment is added, it seems 
that Tx has been broken. However, the system still has a 
special helical or "screw" symmetry: combine the trans- 
lation along X with a rotation, by defining 

Z=Tx{a)Rx{aq), (15) 

where a is the lattice constant and Rx{0) is the rotation 
transformation around axis x. We can see the system is 
invariant under T^, Ty and T^, and these three operators 
all commute. When we turn on the electric field in yz 
plane, this spiral translation symmetry is broken by the 
external field, and we have a charge density modulation 
which is forbidden if the helical symmetry is unbroken., 
and the period is the same with the spiral moment con- 
figuration. However, if the electric field is along x, this 
spiral symmetry is not broken as long as the external 
field is homogenous. We can conclude that, although 
our calculation above was based on an assumption about 
the three energy scales in the problem, this CDW effect 
should be relatively independent of these assumptions, 
because it is controlled by symmetry breaking. 

This absence of accumulation in the absence of an ap- 
plied electric field is strictly correct in the continuum 
limit where the helical transformations are exact symme- 
tries: lattice effects can lead to a very small accumulation 
if the helix vector is commensurate with the periodicity 
of the lattice. In the continuum, charge accumulation is 
always weak as long as the magnitude of the magnetic 
ordering is constant, because of a local symmetry. 

This absence of accumulation in the absence of an ap- 
plied electric field is strictly correct in the continuum 
limit where the helical transformations are exact symme- 
tries: lattice effects can lead to a very small accumulation 
if the helix vector is commensurate with the periodicity 
of the lattice. In the continuum, charge accumulation is 
always weak as long as the magnitude of the magnetic 
ordering is constant, because of a local symmetry. First 
write and (|SJ) in operator form: 

+/jm(f) • i^Pis^pijp). (16) 

Now in this Hamiltonian let the local magnetization m 
vary slowly from point to point. The direction of local 
moments in the whole material can be made uniform by 
performing a position-dependent rotation, written as 

R{r),,m{r)j = m(0),. (17) 

We can keep the interaction H' invariant if we perform 
another position-dependent SU{2) transformation U{f) 



on the spinor ipa, as long as 

UiryS,U{f)^Rir)^,S,. (18) 

Normally, this kind of symmetry is broken by the kinetic 
term because the derivative will be modified to di — > 
di — WdiU after a position dependent transformation. 

However, because we assumed the variation of local di- 
rection is slow, both Rij and U are very slowly varying 
functions, so we can ignore U^diU in the kinetic term. 
This is equivalent to ignoring derivatives of the SU{2) 
transformation. Because the AHE only involves particles 
close to the Fermi surface, as long as the characteristic 
length scale of variation of local moment is much larger 
than 1/kf, the derivatives U^diU do not modify the spec- 
trum close to the Fermi surface and can be neglected. 

Because of this local SU (2) symmetry, the slowly vary- 
ing system is equivalent to the case with homogenous 
magnetization, and without electric field the charges are 
homogenously distributed. However, when the magni- 
tude of local moment varies, the SU (2) gauge symmetry 
is broken as we can never perform any position depen- 
dent rotation to transform this system to the homogenous 
magnetization case. When we add certain external fields 
to the system with constant |m|, the gauge symmetry is 
also broken: local rotation transforms the external field, 
and without protection by symmetry, there will generi- 
cally be an inhomogenous charge distribution. 

In the helical case, an electric field along the helix di- 
rection preserves the symmetry and the charge density 
remains homogenous. It is interesting to look for mag- 
netic moment configurations in which the charges remain 
homogenous for any applied electric field. The constant 
magnetic moment is certainly an example, but there are 
others: to prevent charge accumulation, the Hall current 
must satisfy V • 3 ah = 0. This gives a constraint on the 
moment configuration: 

V • J AH = CT^ffV • (m X E) 

= crAH(E ■ V X m - m • V X E). (19) 

The second term in the bracket of the second line is zero 
because curl of electric field is always zero. Hence there 
is no charge accumulation in any electric field if m has 
zero curl. This requirement is satisfied if m = V<^, so the 
moment configurations with no accumulation are equiva- 
lent to source-free electric field configurations of constant 
magnitude. One example is the monopole configuration, 
in which m(r) = f with the potential 4> = I''!- For this 
monopole background, the charge distribution remains 
homogenous for any applied field. 

We conclude that the coupling between itinerant elec- 
trons and local moment gives rise to an instability in the 
charge distribution in the presence of an applied elec- 
tric field. The simplest steady state of this nonequilib- 
rium problem has an inhomogenous charge density that, 
for a helical magnet, is a charge density wave of the 
same period as the helix. In general, magnetic configura- 
tions with constant magnitude but slowly varying direc- 
tion give rise to charge accumulation only in an electric 
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field, while configurations with variable magnitude cause special configurations where charge homogeneity is pro- 
charge accumulation in equilibrium. However, there are tected for any an applied field. 
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